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I. INTRODUCTION 

The abihty of fabricating artificial materials with so- 
phisticated electromagnetic properties have generated 
large interest recently. By tailoring the electric permit- 
tivity and the magnetic permeability, one can obtain 
arbitrary metamaterial-induced coordinate transforma- 
tions, leading to devices such as perfect lenses and elec- 
tromagnetic cloaks [1, 2, 3]. Unfortunately, the perfor- 
mance of these devices are strongly limited by loss. To 
overcome this limitation it has been suggested to intro- 
duce gain [i[Ei,l3. 

It is crucial to understand fundamentally how active 
media behave electromagnetically. Active media are not 
new; lasers and optical amplifiers are well-established 
technologies. The new aspects are rather the flexibility 
of metamaterials and their applications. For lasers, gain 
saturation plays a vital role when modeling the medium. 
Nevertheless, in applications such as optical amplifiers, 
the medium behaves linearly; thus the linear suscepti- 
bilities and the complex refractive index characterize its 
operation. For potential applications to metamaterials, 
perfect lenses, manipulating the near-fields, cloaks etc., 
linear operation is also desired. Nevertheless, gain sat- 
uration may for example limit the resolution associated 
with the perfect lens, so nonlinearity needs in many cases 
to be taken account of. 

Here we will limit the discussion to linear media. In 
other words, the fields are nowhere allowed to be larger 
than the threshold for gain saturation. Provided there 
are no instabilities, this may be achieved by ensuring 
that the excitations are sufficiently weak. (The fields 
are however considered classical such that the stimulated 
emission dominates the spontaneous emision.) On the 
other hand, if there are instabilities, a linear, frequency- 
domain model is certainly no longer valid. However, if we 
assume that the medium is dark at some time t = Q (for 
example by having the power supply turned off at i < 0) , 
a linear transient analysis is still valid until the fields have 
grown above the gain saturation threshold SgJ. This is 
the standard approach to analyze active, possibly un- 
stable systems in other contexts such as electronics and 
control engineering. In the present work, the transient 
analysis is central for understanding the properties of a 
semi- infinite medium (Section IIIip . 



Recently there has been some confusion about how to 
determine the sign of the refractive index (or the direc- 
tion of the wave vector) in active media d, [^, [13, [HI [H, 
[isl . [31 • The confusion in the field is probably due to 
a nice property of passive media: One does not have to 
invoke causality in its most primitive form to determine 
the sign of the refractive index. One can determine its 
sign straightforwardly at a single frequency by requiring 
that the Poynting vector points away from the excita- 
tion, or by requiring that the wave decays as it propa- 
gates away from the excitation. This makes it natural 
to consider general methods for identifying the sign from 
the (relative) permittivity e(w) and permeability pl{uj) at 
a single frequency lo, even for active media. While such 
methods may give correct answer for a restricted class of 
media, they are necessarily incorrect in general. Indeed, 
there exist both positively and negatively refracting, non- 
magnetic media with e(a;) = 1 — 2ia, < a ^ 1, at a 
single frequency [T5| : Conventional gain media refracts 
positively, while the right-handed negative index medium 
suggested by Chen et. al. [l^, refracts negatively. To 
find the sign of the refractive index, we must go back 
to first principles, using causality in its most primitive 
form: The front velocity of an electromagnetic wave can- 
not travel faster than c, the vacuum velocity of light. (Ac- 
tually, the correct solution can be found without invoking 
causality as an extra principle, as causality is built into 
the Maxwell equations.) Similarly to the classical treat- 
ment of passive media [ITil, this amounts to requiring the 
refractive index to be analytic in some upper half-plane 
of complex frequency, and requiring n— !--)-lasaj— >oo 
[l5| . However, one must watch out for absolute insta- 
bilities; the half-plane may be Imcj > 7 > instead of 
Imo; > 0. 

The purpose of the present article is five-fold: First 
we will consider microscopic causality of active media 
(Section |Tl|. As opposed to passive media, it turns out 
that causality is not sufficient for establishing Kramers- 
Kronig relations. This has to do with possible instabili- 
ties. Second, we will clarify how to resolve the refractive 
index and the direction of wave vector, and pay particu- 
lar attention to the ambiguity mentioned above (Sections 
mil and |V|) : Two media with identical e(cj) and /i(a') at 
LD may respond completely different to monochromatic 
excitations at w. In particular we analyze a slab of thick- 
ness d in 2D. Assuming stability, a nonmagnetic medium 
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with e{w) = \ — 2ia at a single frequency (0 < a <C 1) will 
refract positively when d is small, and negatively when d 
is large. In other words, all media that make the slab sta- 
ble for small d refract positively at the interfaces, while 
all media that make the slab stable for large d refract 
negatively. Third, we discuss and categorize the different 
instabilities that may arise in active media (Section [IV|, 
and relate them to previous literature in plasma physics. 
Fourth, we will prove that for active media, any refractive 
index function in a finite bandwidth can be realized ap- 
proximately, with any precision (Section IVip . Fifth, we 
consider the ultimate limits of active media, and show 
that in general, there are no upper limit of the gain asso- 
ciated with media without absolute instabilities (Section 
IVII|) . Also we show that there are no lower bound for 
the (maximal) gain associated with right-handed, nega- 
tive index media. Permittivity functions of such media 
with arbitrarily low gain are provided, giving directions 
to practical realizations. 



II. CAUSALITY FOR ACTIVE MEDIA 

As mentioned above, we restrict ourselves to linear me- 
dia. Furthermore, the media are assumed isotropic, ho- 
mogeneous, and without spatial dispersion. We also as- 
sume that the medium is time-shift-invariant, that is, if 
an excitation is shifted in time, the response shifts by an 
equal time. Consider the case where the applied electric 
field is due to a external source with charge density Pe (0 
[m . The effect of the source is to induce a charge den- 
sity Pi(i), which together with Pe{t) makes up the total 
charge density p(t) — Pe{t) + Pi(t). Since the medium is 
linear and shift-invariant, we have the following relation 
between the source and the induced charge density: 



Xe{T)pe{t - T)dT. 



(1) 



Here Xe{t) is a real response function, and the lower limit 
in the integral is due to causality; the induced charge 
density cannot precede the source. We will assume that 
the induced charges do not blow up faster than exponen- 
tially for any bounded excitation, which can be formu- 
lated as follows: 



\xe{T)\dT < Cexp(7i), 



(2) 



for some nonnegative constants C and 7. Transform- 
ing ([1]) to the frequency domain, we can write pi{uj) — 
Xe{Lo)pe{oj) and therefore 



p{uj) = [1 + Xe{ljj)]pe{uj), 



(3) 



in an obvious notation. Note that for active media in 
general, the transformation must be performed with the 
Laplace transform, since the induced charge density may 
blow up with time. This amounts to requiring all quan- 
tities to vanish for negative time, and putting Imcj > 7 



in the integrals. For example, 

P^{lo) = I pe{t) exp(iujt)dt. 



(4) 



Eq. ^ together with the Maxwell equations eo V • E = p 
and V • D = pe, and the constitutive relation D = eeoE, 
mean that 
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1 + Xe{Lj) ' 



(5) 



— ^ = 1 + / Xe{t) exp{iujt)dt. (6) 
Jo 



Since Xe{t) is real we also have the symmetry 
e(-tj*) = e*(w). 



(7) 



We note that it is not the permittivity itself, but its re- 
ciprocal, that is causal in the sense of being the transform 
of a response function that vanishes for t < [l^ . 

At very high frequencies, the electrons behave essen- 
tially as if they were free, yielding the asymptotic form 



m 



00. 



(8) 



While causahty (O, symmetry ([7]), and the asymptotic 
form ([5]) are valid in general, Kramers-Kronig relations 
cannot be established unless more information on the me- 
dia is known. In other words, even though media that 
satisfy the ususal Kramers-Kronig relations are causal, 
not all causal media satisfy the Kramers-Kronig rela- 
tions. This is particularly true for active media, which 
may exhibit certain instabilities. For example, the re- 
sponse function Xe(t) may increase exponentially, which 
leads to a singularity of l/e(a;) in the upper half-plane 
Imw > 0. Or the transformed response function Xe{i^) 
may be equal to —1 somewhere in the upper half-plane, 
giving a singularity of e{uj) there. 

In light of Titchmarsh' theorem [l^ [23|, if no such 
singularities are present, not even at the real frequency 
axis, (O and ([5]) give the Kramers-Kronig relations 



Ime(w) = — V 



Ree(w') - 1 



duj\ 



2 f°° Ime 
Ree{u;)-1 = -V — 



Ime(a;')cj' 



duj\ 



(9b) 



where V denotes the Cauchy principal value [331 ■ If 
e(w) has singularities at the real axis, the Kramers- 
Kronig relations must be modified. For example, if the 
medium is conducting at zero frequency, e(a;) is singular 
at cij = 0. Then, the Kramers-Kronig relations are re- 
tained if we subtract the singularity, i.e., make the sub- 
stitution e(u;) — > e(w) — ia/uj in ([5]), where > is the 
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zero frequency conductivity [2]| . If e(u;) have singulari- 
ties in the upper half-plane, e(a;) loses its meaning at real 
frequencies and the Kramers-Kronig relations cannot be 
expressed along the real axis. The observation frequency 
uj and the integral must instead be taken along a line 
above the singularities. We will discuss media with in- 
stabilities further in Section [TVl 

We can treat the magnetic permability in a similar 
fashion, yielding 



1 



Xm{t) ex.p{iu!t)dt, 



(10) 



instead of Here Xm{t) is a real response function. 
With a similar asymptotic form as ([8]), and in the ab- 
sence of singularities in the upper half-plane, we may also 
obtain Kramers-Kronig relations for the magnetic per- 
meability. Using a Kramers-Kronig relation analogously 
to (j9bp to calculate Re/i(0) from Im/i([j), we straight- 
forwardly find that Re/i(0) > 1 whenever Im/i(tj) > 0. 
This is not always true |22] , as demonstrated by the exis- 
tence of passive, diamagnetic media. In other words, the 
Kramers-Kronig relations for the magnetic permeability 
must be treated with care. 

For the remaining parts of this paper, we will only as- 
sume causality in the sense © and (fTU]) . meaning that 
we allow for instabilities. However, to retain the mean- 
ing of complex frequencies, we exclude superexponential 
instabilities; in other words ((2) and a similar bound for 
Xm{t) are required d^]. Similarly to ([5]) we will assume 



e(a;),/i(a;) ^ 1 as a; ^ +oo. 



(11) 



whenever real frequencies are meaningful, and Eqs. ^ 
and (jlOp imply 



e(w),/i(w) 



1 as Imw — > +00. 



(12) 



III. 



DIRECTION OF THE WAVE VECTOR AND 
SIGN OF THE REFRACTIVE INDEX 



To identify the direction of the wave vector and the 
sign of the refractive index in active media, one can use 
the following approach : By a transient Laplace trans- 
form analysis of a slab of thickness d, one may compute 
the fields for time t < d/ c, where c is the vacuum velocity 
of light. Then, by causality, the fields have not felt the 
presence of the far end; thus they must be identical to 
those of a semi-infinite medium. By subsequently taking 
the limit d — s- oo, one can extract the wave vector or re- 
fractive index. Below we will obtain the same results as 
in ^5*1 more directly, by taking the limit d — > cxd imme- 
diately in the region of convergence of the Laplace trans- 
formed fields. Consider a slab with permittivity e{uj) and 
permeability /x(a;), surrounded by vacuum, see Fig. [1] 
Since the medium may be active, the fields may blow up 
with time. Consequently, steady-state frequency-domain 
fields do not necessarily exist. The natural remedy is to 




FIG. 1: (Color online) Slab of thickness d. 



consider Laplace transformed fields instead of nonexist- 
ing Fourier transforms, assuming the field is equal to zero 
for negative time: 



E{uj) = I £{t)exp{iujt)dt. 
Jo 



(13) 



Here £{t) denotes the real, physical, time-domain elec- 
tric field. The Laplace transform exists in a region 
of convergence Imw > 7, where 7 is a sufficiently large, 
positive parameter, such that the exponential factor en- 
sures convergence of the integral. The inverse transform 
is given by 

£{t) = —j E{uj)eyip{-iujt)duj. (14) 

Let the Laplace transformed, incident wave at z = 0^ 
be E^{u)) = Eq {lo) ey^Y'i'ikxX iky-y), where Eq{uj) is 
independent on the spatial coordinates, and and ky are 
the transversal, spatial frequencies of the source. Solving 
Maxwell's equations in the transform domain, we obtain 
the field 



E{lu)/E+{lu) 

exY>{ikzz) + Rexp{—ikzz) 

exp{ik'^z) + S~ exp(— i/c^z) 
T exj>{ikz{z — d)) 



(15) 

for z — 0^ , 
for < z < d, 
for z > d. 



where the total reflection coefficient R, field amplitudes 
in the slab [S'^ and S^), and transmission coefficient T 
are found to be 



S+ 



s- 



(p^kl-ki')[l-exp{2tk'J)] 



(/ifcz 



k'J^ - (pkz - k'^f exp{2ik'J) 
2^ikz{iJ.kz + k'^) 



{^ikz + k'^Y - {^jikz - fc^)2 exp{2ik'J) ' 
2^fc^(/ifc^ - fc^) 

' ipkz - k'^y - i^ikz + k'zY exp{-2ik',d) ' 
4^kzk'z exp{ik'zd) 

{fikz + k'^y - {pkz - k'^Y exp{2ik'^d) ' 



(16a) 
(16b) 
(16c) 
(16d) 



Here, k^ = w^/c^ — k^ — ky and fc!,^ = tixuY jc? — 
The expressions above apply to TE polarization; for TM 
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polarization, the expressions are valid provided e and fi 
are interchanged. Note that the fields are invariant if 
k'^ —k'^, so the choice of the sign of k'^ does not matter. 

For active media, i?, S, and T may contain singularities 
in the upper half-plane. Thus the fields should only be 
evaluated in their region of convergence Imo; > 7, where 
the line w = 17, 7 > is located above all singularities 
[23I I . Note that singularities in the upper half- plane have 
nothing to do with noncausality; they rather imply insta- 
bilities In practice, diverging fields imply that the 
slab starts lasing (taking gain saturation into account). 

To obtain the fields associated with a semi-infinite 
medium, we may take the limit d ^ 00. From p2|l . it 
is clear that k'^ k, -iiuj/c for sufficiently large Imw. Since 
the sign of fc^ is arbitrary in (fT6)) . we take k'^ k, ~\-uj/c 
there. With this choice, Im/c^ is positive for sufficiently 
large Ivauj. It is not difficult to realize that the singular- 
ities of S do not move towards oj = ioo as the thickness 
d ^ 00; therefore we may take the limit d ^ 00 to find 
the field in a region Imcj > 7. The resulting expression 
is given by 



the simple formula 



S = 



^kz 



— exp(ifc^z). 



(17) 



We have argued for (fT7)) for Imw > 7. The real, time- 
domain field is obtained by inverse transforming SE^ [uj) 
along the Bromwich path uj = ij. However, by analytic 
continuation and Cauchy's integral theorem, we may 
move the integration path down towards the first non- 
analytic point of (jl7[) ; in other words we may now set 7 = 
max{0, Im (zero of fikz + k'^), Im (nonanalytic point of k'^)} 
In many cases of interest, we obtain 7 = 0, which means 
that (|17p and k'^ may be interpreted at real frequencies. 
In these cases, k'^^Lo) is the analytic continuation of 

the branch of y^e/iw^ _ _ ^2 tends to +iLj/c as 

l^l — > 00. Thus if possible, the sign of k'^ must be chosen 
such that k'^ is analytic in the upper half-plane, and 
such that k'^ — s- +lo/c as \uj\ ^ 00. If this is not possible 
due to poles or odd-order zero(s) of e/iw^ — k^ — ky, 
there are instabilities. Then k'^ does not have any 
physical significance for real frequencies. The nature of 
instabilities associated with poles or odd-order zeros of 
e/jLUj'^ — k. 

Instabilities associated with zeros of fik 
to the boundary of the medium, and may be eliminated 
by choosing another surrounding material [15 1 



ky will be discussed in the next section. 

fc' are related 



n{u;) = v^7MR^exp[i(^,(^) + <Pp(^))/2)], (18) 

where (pt{Lo) -f (Pfj,{Lo) is the complex argument of 
e{uj)fi{uj), unwrapped such that it is continuous for lo ^ 
and such that it tends to as w ^ ±cxd. If e(uj)fj,{uj) does 
have poles and/or odd-order zeros, there will be instabili- 
ties, and n{uj) does not have meaning at real frequencies; 
only above the nonanalytic points. Note however, that 
relativistic causality never is violated [isj : n(w) is always 
analytic in some upper half-plane Imw > 7 > 0. 

Due to the phase unwrapping procedure, the sign of 
n{uj) at a particular, real frequency is determined from 
the global properties of the functions e(cj) and fi{uj) 
[Tsl fig . UJ, [2^ . As demonstrated in Section |Vl two ma- 
terials with the same permittivity and permeability at a 
particular frequency, may have opposite sign of the re- 
fractive index. Therefore, any method that identifies the 
sign of n{u) from e(w) and fJ-{uj) at a single frequency 
u) must he incorrect in general. Consequently, although 
the methods for identifying the sign of n(uj) in Refs. 

d, i, [M [H E [S [3 may give the correct result for 
certain active materials, they are incorrect in general. 

A critical point which cannot be overemphasized, is 
that the complex frequency-domain fields are not the 
physical fields themselves; they are suitable transforms 
of the real, physical, time-domain fields (Fourier trans- 
forms for passive systems, and Laplace transforms for ac- 
tive systems). Laplace transformed functions exist only 
in their region of convergence Imw > 7; thus the trans- 
formed fields must not be interpreted elsewhere. For ex- 
ample, it is ternpting to take d ^ 00 while letting oj 
be real in (jl6p [12l |. However, for media that give rise 
to singularities of Ijl6p in the upper half-plane, this is 
clearly not a correct procedure and may in fact result 
in incorrect sign for the refractive index. This happens 

e. g. for a conventional gain medium such as the inverted 
Lorentzian medium. 

Finally, we should keep in mind that any physical 
medium is limited by gain saturation. When there are no 
instabilities, the linear model is valid provided the excita- 
tions arc sufficiently weak. When there are instabilities, 
the time-domain fields are only valid before they reach 
the threshold for gain saturation. Provided this happens 
sufficiently late, transform-domain concepts such as the 
refractive index and the wave vector have meaning. 



Taking kx 



0, we have fc^ = n{u})uj/c, where 
n{u!) — ±-\/ e{uj)fi{Lu). Thus, provided e{Lu)fi{u!) does not 
have any poles or odd-order zeros in the upper half-plane, 
n(w) is identified as the analytic function in the upper 
half- plane that tends to -1-1 as w ^ 00. This result was 
used in Ref. [l6l | to determine the refractive index in 
certain active, nonmagnetic, negatively refracting media. 
When n(Lu) can be identified as an analytic function in 
the upper half-plane, and e(tj)/i(tj) is continuous at real 
frequencies (except possibly at w = 0), n{Lo) must be 
continuous at real w 7^ 0. Thus we can identify n(u;) by 



IV. INSTABILITIES OF CAUSAL MEDIA 

An instability means that the electromagnetic fields 
blow up with time. There are three different classes of 
instabilities [1^ [13, HI] . First we let the medium be in- 
finite. An absolute instability means that the fields at 
a fixed point in space blow up with time. A convective 
instability means that the fields blow up with time, but 
at a fixed point, they do not. This means that the in- 
stability is "convected away" . Absolute and convective 
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instabilities are indicated in Fig. [2] A third category, 
global instabilities, arises when the medium is bounded. 
For example, the fields may diverge due to amplified, 
multiple reflections. 




convective instability z 



FIG. 2; Absolute and convective instabilities. 

An example of a convective instability is that of con- 
ventional gain media, such as the Erbium-doped fiber 
amplifier. When such a gain medium is put into a res- 
onator configuration, a global instability may occur. Due 
to gain saturation (which is a nonlinear process), this 
leads to lasing rather than infinite fields. 

Absolute and convective instabilities are features of the 
medium itself, not the total system including boundaries. 
Absolute instabilities arise when the refractive index 
cannot be identified as an analytic function in the upper 
half-plane Imw > 0. This happens if e{u!)fi{Lu) contains 
poles or odd-order zeros there. 

An apparent absolute instability arises at oblique in- 
cidence to a semi-infinite medium (25j . Even if n{uj) is 
analytic in the upper half-plane, k'^ at oblique incidence 
may have branch points there. This can be interpreted as 
follows: Any causal excitation contains necessarily an in- 
finite band of frequencies. Some of the frequencies make 
e/iw^ — k'^ — ky zero; this corresponds to waves propa- 
gating perpendicularly to the z-axis. Such waves propa- 
gates an infinite distance, and therefore pick up an infi- 
nite amount of gain, before arriving the plane z — const. 
This is called a virtual absolute instability, as it appears 
mathematically as an absolute instability, while its phys- 
ical interpretation is similarly to convective instabilities. 
Virtual absolute instabilities can be eliminated by lim- 
iting the extent of the active medium in the transverse 
direction. 



V. STABILITY AND REFRACTION IN 2D 

We will now limit the discussion to absolutely stable 
media for which e{uj) and fi{uj) do not have singularities 



or zeros in the upper half-plane. (We do not, however, 
exclude the possibilities of convective or global instabili- 
ties.) 

The uniqueness theorem for solutions to Maxwell's 
equations in the frequency domain [29] assumes that the 
medium is lossy. Therefore, it is not surprising that 
for active media, e and at a single frequency u do 
not determine whether the material refracts positi vely 
or negatively at w. Indeed, as pointed out earlier p^ . 
there are both positively and negatively refracting me- 
dia with identical electromagnetic parameters at a single 
frequency. 

For a semi-infinite material the sign of k'^ determines 
the reflection, refraction, and propagation of an incom- 
ing wave. In a slab the sign does not matter as the 
fields ([TB)) are invariant under k'^ ~k'z- This seems 
to be a contradiction: If a beam is incident to the slab 
at some angle of incidence 9, the first-order transmit- 
ted beam must exit the slab either above or below the 
entrance point, see Fig. [31 A natural question arises: 
For a slab of thickness d, and e = 1 — 2ia, < a <C 1 
and /i = I at the excitation frequency, will the beam 
refract positively or negatively at the interfaces? To in- 
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FIG. 3: (Color online) Schematic view of a ray with oblique 
incidence to a nonmagnetic slab with e = 1 — 2ia. 

vestigate this problem, we consider active materials with 
e = 1 — 0.030* and = 1 at a normalized frequency 
uj — 1. We also assume that for any d, there is such a 
material that makes the slab electromagnetically stable 
(no global instabilities). This assumption will be justified 
below. Let the frequency-domain excitation be a beam at 
oblique incidence. The beam comprises a superposition 
of plane waves, whose amplitudes A(kx) are Gaussian 
distributed: A{k^)_^ l/Vina"^ exp[-{k^ - h)^/{2a^)] 
where a = 0.032, kx = 0.50, and < fc^: < I, normal- 
izing the vacuum light velocity [c — 1). To determine 
where the beam exits the slab, we consider the function 
A{kx)T{d,kx), where T = T{d,kx) is the transmission 
coefficient (I16dp . By evaluating the inverse Fourier trans- 
form in kx, we get the field E{d, x) at the right-hand in- 
terface. Thus we can determine the correct solution (pos- 
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itively or negatively refracted beam) for different d. The 
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FIG. 4: (Color online) The field amplitude |£;(d,a;)| (arbi- 
trary units) as a function of slab thickness d and transversal 
coordinate x. 

field amplitude \E{d,x)\ is plotted in Fig. HI The beam 
undergoes multiple reflections in the slab; thus it exits the 
slab at several different locations along the x-axis. We 
observe that the beam refracts positively for d ^ 250, 
and negatively for d 290. For 250 ^ d ^ 290 the re- 
sults indicate that there are no media with e = 1 — 0.03i 
that make the slab electromagnetically stable. 

The results certainly do not imply that for a given 
medium, the beam refracts positively or negatively de- 
pendent on the slab thickness. They rather imply that 
the media that make the slab stable for small thicknesses 
refract positively, while the media that make the slab 
stable for large thicknesses, refract negatively. This is 
indicated in Fig. [51 The first category of media sup- 
ports an amplifying, forward propagating wave, that is 
refracted positively. For this category, when d ^ 250 the 
round-trip gain becomes larger than the Fresnel losses, 
and the slab is unstable (global instability). The second 
category supports a backward propagating wave that is 
refracted negatively, and amplified as it propagates to- 
wards the left interface ^30]. For this category, a global 
instability arises when d ^ 290. When the slab thickness 
becomes large, the amplitude of the backward wave is 
small at the right-hand interface. This explains the de- 
caying field amplitude as d is increased above 290 in Fig. 

H 



unstable 


^^^^■^ n = —l-\-ia 







"limi "limj Slab thickness d 

FIG. 5: (Color online) Stability of nonmagnetic Id slabs with 
e = 1 - 2ia. 

We now concretize the two categories of media by 



examples. Consider two nonmagnetic materials with 
ei(w) = 1 + Xgiuj) and 62(0;) = [1 -I- Xr(w)]^, respectively, 
where Xr,gi^) are Lorentzians in the form 



F 

Lo'^ g - Up- - iL0Tr,g 



(19) 



With the parameters Fg = -0.00016, Tg = 0.005, and 
ujg = 1 for material 1, and F,- ~ 3.7, F^ = 0.005, 
and LUr = 0.6 for material 2, these two materials have 
£1,2(^3) = 1 - 0.030i. With the help of (HI]), we find 
the refractive indices ni{u!g) w 1 — 0.015« and n2{ujg) ~ 
— 1 + 0.015*. These materials are both right-handed. 
They do not support absolute instabilities as ei^2i'-^) have 
no poles or zeros in the upper half-plane. Clearly, both 
of them show convective instabilities. 

We should now investigate the stability as a function 
of d, when the slab is made of either of the two example 
materials. First we limit ourselves to a ID slab. Active 
materials in ID may for example be realized in the form 
of a transmission- line model with lumped circuit elements 

fc„ = in (Ull, 



[30|. For ID propagation, we take = h 
and determine numerically when there are no poles in the 
upper half-plane. We find that for d < dumi ~ 300 the 
material 1 slab is stable, while for d > duma ~ 320 the 
material 2 slab is stable. 

In the 2D case stability requires that there are no poles 
in the upper half-plane of (|16|) for any k^- For material 2, 
it turns out that the stability limit remains dum^ « 320. 
On the other hand, for material 1 the stability limit turns 
out to be zero (dumi — 0). By reducing the slab thick- 
ness, (|16p may not have poles for a fixed angle of inci- 
dence, but there will always be poles for a sufficiently 
large angle of incidence. This is due to the fact that as 
the angle of incidence increases towards 7r/2, the wave 
propagates longer between two refiections. Thus, during 
one round-trip, the wave picks up more gain compared to 
the losses associated with the two reflections. By limiting 
the transversal dimensions of the slab, and limiting the 
reflections at the upper and lower boundaries (e.g. us- 
ing anti- reflection coatings or absorbing layers), this 2D 
instability can be eliminated. 

As we have seen, the fields of a slab are invariant to 
the sign of the refractive index. This is also the case for 
several other propagation problems such as plane wave 
scattering on a sphere or wave propagation in waveguides 
[sH . However, the fact that the sign is invariant does 
not necessarily mean that it is irrelevant. The refrac- 
tive index indicates the direction of the wave ( "forward" 
or "backward") before multiple reflections dominate the 
picture; it gives you information on whether the wave will 
refract positively or negatively at a boundary etc. Thus 
one can determine several properties of a given system by 
evaluating pSl) . without solving Maxwell's equations in 
detail. Nevertheless, for active materials it is important 
to be aware of possible instabilities, which may make 
frequency-domain concepts useless for real frequencies. 
The only way to fully understand whether a given shape 
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of an active material is stable or not, is to solve Maxwell's 
equations with appropriate boundary conditions. 

VI. LIMITATIONS FROM CAUSALITY 

Passive media, i.e., media in thermodynamic equilib- 
rium in the absence of the variable field 21], satisfy 

Ime(w) > and Im/i(w) > for w > (20) 

in addition to This leads to fundamental limi- 

tations for left-handed media; for example, there is a 
lower bound for the dispersion associated with transpar- 
ent, left-handed media, and there is a lower bound for 
the loss associated with nondispersive, left-handed me- 
dia [HI,!!!!!!. 

Active media do not have the limitation (|20|) . i.e., 
Ime(a;) and Im I-J,{lu) may take any value. It is interesting 
to see if there are still any limitations resulting from sym- 
metry ^ and causality in the sense ([9]). It turns out that, 
in fact, on a finite bandwidth, there are no fundamental 
limitations. For example, e(w) = /x(aj) = — 1 (lossless 
Icft-handedness) can be approached on a nonzero, finite 
bandwidth. This result follows directly from a standard 
result of the theory of Hardy spaces: A function e{uj) 
satisfying ((Tj) and ([9]) can approximate any square inte- 
grable function /(lu) on a finite bandwidth 41]. The ap- 
proximation may be achieved with any precision; in the 
Krem-Nudel'man case 34] at the expense of the norm 
of €{uj) — 1 outside the bandwidth of interest. Of course, 
a large norm outside the relevant bandwidth may im- 
ply difficulties of realization. Nevertheless, contrarily to 
claims elsewhere [35], we note that causality does not 
prohibit e.g. e{uj) = fJ.{uj) ~ —I with any precision, even 
in a finite bandwidth. 

The possibility of approximating any desired behav- 
ior in a limited bandwidth may seem useless unless 
the resulting medium is free from absolute instabilities. 
With a straightforward approximation using e.g. Krcin- 
Nudel'man, the resulting permittivity may not be zero- 
free in the upper half-plane. A possible remedy is to 
approximate log/(ti;) with a function g{uj) using Krem- 
Nudel'man, and setting e{uj) = exp[g{u!)]. With this ap- 
proach, it is necessary to assume that log/(ix') is square 
intcgrablc; in other words, /(w) is not allowed to be zero 
on an interval of nonzero measure. This procedure clearly 
gives an analytic and zero-free function e{io) in the upper 
half-plane. Square integrability of e(a;) — l (and Kramers- 
Kronig relations) are satisfied provided log/(aj) is Lips- 
chitz continuous and approaches zero at the end-points 
of the interval 36] . (The latter condition is not a further 
constraint as it may be fuUfilled by extending the original 
interval.) 

As an example of gain compensation of the losses asso- 
ciated with a left-handed resonance, consider the causal 
medium e(aj) — 1 -I- Xr(w) -I- Xg(w) and //(w) = l + Xri^), 
where Xr,g{^) are the Lorentzians defined in (I19p . Tak- 
ing ujr = I, uig = 1.058, Fr = 0.25, Fg = -0.0034, 



= 0.005, and Tg = 0.02, we find that n(w) = -1 
and dlmn^(uj)/duj = for uj — 1.06. This medium has 




Re n 



0.8 1 1.2 1.4 1.6 

Frequency uj 

FIG. 6: (Color online) Refractive index for the left-handed 
medium with gain compensation. At = 1.06 we have lossless 
negative refraction, n = — 1 and dlmn'^ {uj)/duj — 0. 

Ime(a;) < in the bandwidth [1.05, 1.07]; thus it is net 
active there. This does not imply that the system is un- 
stable unless the medium is infinite or put in a resonator 
configuration; examples of such stable systems include 
fiber optic amplifiers. As argued in Section Hvl the ab- 
sence of absolute instabilities of this medium is guaran- 
teed by the fact that e{uj) and n{uj) do not contain zeros 
or poles in the upper half-plane. We note that the gain 
compensation due to Xg(<^) has completely removed the 
loss which would have been present in the absence of this 
gain. 

Lossless left- handedness in a finite bandwidth (0^1,072) 
(with < LUi < UJ2 < 00) may in fact be obtained with a 
rather general class of susceptibility functions, even if we 
restrict ourselves to passive media. Pick any real, square 
integrable function /(w), with f{uj) = in {uji,ijJ2)- Set 
Ime(a;) = Im/i(a;) = fiuj). Thus, by tailoring /(w) this 
medium is defined lossless at and near w, and may be de- 
fined lossy otherwise. Identifying Ree(a;) = Re/i(a;) with 
the Kramers-Kronig relation (|9bp , the resulting medium 
is causal by definition. From 

Ree(tj) = Re^(w) (21) 
_ ^ 2uj'f{uj')dLj' r°° 2uo'f{uo')duj' 

Jo w2 - J^^ _ ^2 

we observe that provided f{uj) is positive and sufficiently 
large below uji, the medium will be left-handed. In fact, 
one can in principle construct causal, passive, negative 
index media with arbitrarily low (maximal) loss (see Ap- 
pendix A). 

VII. BOUNDS FOR THE GAIN 

The permittivity of an inverted Lorentzian medium, 
ei{io) = 1 -|- Xg(w), has a zero in the upper half-plane 
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when Fg < —1. Thus a sufficiently strongly pumped 
medium seems have an absolute instability. It is natural 
to ask whether there is a general, upper bound on the 
gain for media without absolute instabilities. It turns 
out that it is not so. In fact, the second material in the 
previous section, £2(1^-') = [1 + Xr{^)]'^ and /i2(i-^) = 1, 
does not have absolute instabilities no matter how large 
Fr and therefore the maximum gain maXij>o[— Ime2(cij)] 
are. This is realized from the fact that 1 + Xr{^) is a valid 
permittivity function of a passive medium, and does not 
have zeros in the upper half-plane. 

Similarly, it is interesting to investigate if there is a 
least maximum gain to obtain a nonmagnetic (right- 
handed) negative index medium. A problem with the 
nonmagnetic negative index media that have been sug- 
gested so far [1^ [s^l, is that they require very high 
gain in some spectral areas adjacent to the working fre- 
quency. This may imply difficulties for realizations, po- 
tential instability problems due to imperfections etc. We 
will now prove that it is indeed possible to reduce this 
gain without destroying the right-handed, negative in- 
dex behavior. In fact, there exist nonmagnetic, causal 
media with arbitrarily low maximum gain, which re- 
fract negatively. To see this, we employ the result of 
Appendix \^ which states that there are causal, pas- 
sive, left-handed media with arbitrarily low maximum 
loss. Examples of such passive media are those given 
by e{uj) = nioj) = 1-1- u(a;) -t- iv{uj), where the required 
properties of the susceptibility u{ijj) + iv(uj) are discussed 
in the Appendix. Letting our active medium instead 
have = [1 -I- u{uj) + iv(ijj)Y and /i(w) — 1, where 
uiuj) + iv{(jj) is the susceptibility of the passive medium 
in the Appendix, n{uj) must be identical to the refrac- 
tive index of that passive medium. From the fact that 
the passive medium can have arbitrarily low maximum 
loss, we realize that the maximum gain of the present 
medium can be arbitrarily low. Thus, we have obtained a 
right-handed, negative index medium with arbitrarily low 
maximum gain. In Fig. [7] the refractive index of such a 
medium is indicated, giving a rough guide to the realiza- 
tion of such media: The permittivity may e.g. be taken 
to be the square of a superposition of narrow Lorentzians. 
While a superposition of Lorentzians might be easy to re- 
alize approximately, it is not necessarily straightforward 
to achieve the square operation. 

For such constructions from passive medium counter- 
parts, (jA3p is valid for any refractive index functions that 
satisfy the Kramers-Kronig relations in the usual form 
([^. In other words, negative refraction requires either 

(i) large imaginary part below the working frequency, or 

(ii) exponentially steep variation immediately below the 
working frequency, or (iii) singularities at real frequen- 
cies. 

For arbitrary active media that satisfy the Kramers- 
Kronig refactions in the usual form, inequality (|A3p re- 
mains valid if there is no gain above the observation fre- 
quency. If there is gain above the observation frequency 




Frequency uj 

FIG. 7: (Color online) Refractive index of a right-handed neg- 
ative index material with low maximum gain. The required 
gain to obtain Ken{uj) = —1 at some frequency can be made 
arbitrarily low by letting the edge of Imn(w) be sufficiently 
steep. The refractive index is n(aj) = 1-1- u{u!) + iv{u;), where 
u{uj) + iv{Lj) is the superposition (|A2|I of Lorentzians. The 
parameters used are loi — Alo = 1 and F = 0.001. 



wi, it becomes 



u{lui) > — 



2Wn 



In- 



(22) 



gam. 



is the maximum frequency with 
l/Aw is the maximum steepness of w(w), and Umax = 
max |Imn([.i;)|. The proof of (|22p goes similarly to that of 
(|A3|) in the Appendix. This shows that there is a trade- 
off between the maximum gain or loss, and the steepness. 

There is however a possibility that negative refraction 
can be obtained with a nonmagnetic medium with large 
maximum loss (in accordance with ^22\i ) but small max- 
imum gain. By letting e(tj) = 1-1- Xr{^) + Xg{^)i where 
Tr = Tg is fixed, it is indeed possible to obtain a refrac- 
tive index with Re n < in some spectral region, while 



\Fg\ is arbitrarily low. However, in the limit Fg 



0, 



the figure of merit — Ren/Imn tends to zero; thus these 
media may not be very useful. 



VIII. CONCLUSIONS 

Several questions related to active materials have been 
addressed. We have shown that Kramers-Kronig rela- 
tions (for real frequencies) cannot always be established 
for causal, active media, due to the possibility of absolute 
instabilities. For active media that satify the Kramers- 
Kronig relations, causality imposes no fundamental limits 
to the dispersion and loss. 

Furthermore, we have considered the direction of the 
wave vector, and emphasized that if possible, the wave 
vector and the refractive index must be chosen analytic in 
the upper half-plane of complex frequency, and such that 
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n — > +1, kz — > +ll}/c as w — > oo. If the wave vector or re- 
fractive index have branch points or singularities in the 
upper half plane, the respective functions do not have 
meaning for real frequencies. Such nonanalytic points 
mean absolute instabilities, which arc fundamentally dif- 
ferent to the convective instabilities of conventional gain 
media, and global instabilities associated with bounded 
systems. 

Materials with e = 1 — 2ia (0 < a ^ 1) and ^ — 1 at 
a single frequency refract either positively or negatively. 
The materials that make a slab stable for small thick- 
nesses, refract positively, while the materials that make 
the slab stable for large thicknesses, refract negatively. 
For a fixed thickness, at most one of the two categories 
of media makes the slab stable. 

Finally we have argued that there are absolutely stable 
media with arbitrarily large gain. Moreover, we have 
proved that there are nonmagnetic negative index media 
with arbitrarily low maximum gain, giving directions to 
realizations. The proof is based on the fact that there 
are passive, left-handed media with arbitrarily low loss. 



APPENDIX A: PASSIVE LEFT-HANDED MEDIA 
WITH ARBITRARILY LOW MAXIMUM LOSS 

Here we will show that there exist causal, passive, left- 
handed media with arbitrarily low maximum loss. The 
proof is by construction. Let the permittivity be written 
e{uj) = 1 + u{uj) + iv{uj), where 1 -|- u(u!) and v{ijj) are 
the real and imaginary parts, respectively. Passitivity 
means that v{uj) > for a; > 0. Since v{uj) — can 
be approached, we allow ourselves to put v{uj) = as 
well; adding a small, slowly varying function to v(oj) does 
not alter the argument below. Let v{io) > vq > for 
luq < Lo < UJl — Alo, and v{uj) — for lu > uJi, see Fig. |8l 
Estimating u{uji) by the Kramers-Kronig relation (j9b[) . 




FIG. 8: The function v{u}) has support below cji. For ujo < 
UJ < UJl — Auj, we assume v{oj) > vo- 



we find 

w(wi) 



< 







v{uj)uj(hj 


71" 


'o 


ujf — Uj'^ 


Vq 




-^"^ 2ujdw 


TT 




ujI — Uj"^ 


VQ 
TT 


L 


UJI-UJQ 

— I hi 

Aw 



2wi 



uJq + UJl 



(Al) 



Now we can choose any small vq and require the maxi- 
mum of v{uj) to be less than, say 2vq. By having a suffi- 
ciently narrow transition region Aw, u(wi) gets less than 
— 1. Setting /i(w) — e(w) = 1 + u{lj) + to(w) completes 
the proof. 

Note that the required function can be approached by 
superpositions of several, narrow Lorentzians, with reso- 
nance frequencies equally spaced in the interval (ujq, uji — 
Aw). For example, in the limit of continuous varying res- 
onance frequencies from wq = to wi — Aw, we obtain 



u{uj) + iv{uj) CC I ^2 2 



WqC?wo 



Wq — w^ — iwF 



(A2) 



= UJl — Auj — s arctan 



UJl — Auj 



where s = i\fu? 



..liviv^ u — I. y I „wF). The required permittivity and 
permeability are obtained by choosing a sufficiently small 
F. 



The steep edge in the transition band may imply that 
the medium is difficult to realize. Given a maximum 
loss i^max (i'('^) ^ Wmax for all w). One can prove that 
such steep edges is the only way to obtain left-handedness 
for media that satisfy Kramers-Kronig relations in the 
usual form ([5]). Indeed, for any square integrable function 
v[uj) > with limited steepness {\dv{uj)/duj\ < 1/Aw for 
some Aw), 



u(wi) > 



In- 



2wi 



,Aw 



(A3) 



Here wi is an observation frequency. The inequality (|A3[) 
is found by a similar argument as that of (|A1[) , consider- 
ing the fact that the least possible w(wi) is obtained when 
f (w) = Wniax for < w < wi — (5w, and f (w) decreases 
linearly to zero above wi — (5w. Here Su is a positive 
parameter. Letting 1 + u(w) + iv{uj) be equal to e(w), 
n{uj), or n(w), the bound (|A3P applies to the permittiv- 
ity, permeability, and refractive index of all passive media 
that satisfies the Kramers-Kronig relations in the con- 
ventional form For media with singularities of e(w) 
or fj,{uj) at real frequencies (e.g. an ideal plasma), (|A3|) 
does not apply. We can conclude that left-handedness 
implies either (i) large loss below the working frequency, 
or (ii) exponentially steep variation immediately below 
the working frequency, or (iii) singularities at real fre- 
quencies. The trade-off between requirements (i) and (ii) 
is quantified by (|A3p . 
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